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Scatter  Channels  and  Target  Modeling  For  Active  Signal  Extraction  in 
Underwater  Acoustics:  II.  (Summary  and  Synthesis  of  I,  [1]) 

by  David  Middleton* 

August  1991 

SECTION  1.  INTRODUCTION 

The  main  purpose  of  the  document  is  to  provide  a  concise  review  and  synthesis  of  the 
detailed  analytic  models  described  in  current  work  of  the  author  [1]  directed  toward  the  active 
acoustic  detection  of  underwater  targets  in  the  ocean.  The  aims  here  are  primarily  descriptive  and 
summarial,  with  the  emphasis  on  motivation  and  interpretation  rather  than  on  the  analytic  details 
which  are  developed  in  [1]  itself  and  which  follow  from  earlier  studies  [2].  Of  course,  an  analytic 
structure  of  sufficient  sophistication  is’ required  in  technical  applications  and  for  the  quantitative 
guidance  of  experiment  In  our  review  we  shall  for  the  most  part  provide  the  main  results  from 
[1]  and  show  how  they  are  important  in  obtaining  effective  detection  of  underwater  acoustic 
targets.  Here,  we  are  concerned  principally  with  (i),  weak  narrowband  signals  in  the  far-field,  and 
(ii)  dominant  signal-dependent  noise  or  reverberation,  since  these  are  usually  the  critical  limiting 
factors  in  performance. 

The  principal  model  components  may  be  identified  typically  as  follows: 

I.  The  generalized  channel  model,  expressed  in  terms  of  (linear)  spatio-temporal  filters  and 
consisting  of  both  scatter  (reverberation)  and  ambient  noise  mechanisms  and  including  the  array 
coupling  at  source  and  receiver; 

n.  Target  models,  which  should  be  analytically  simple  but  still  retain  the  principal  features  of 
typical  targets,  e.g.,  cross-section,  specular  and  diffuse  scattering  mechanisms,  spatial 
distribution,  and  Doppler; 

HI.  Statistical-physical  descriptions  of  the  ambient  and  scatter  components  of  the  channel,  and 
the  target,  e.g.,  moments,  covariances,  and  pdfs,  as  needed  for  the  signal  processing  required 
in  IV  ff;  and 

IV.  Detection  algorithms  and  performance,  whose  parameters  are  determined  by  the  model 
elements  of  I-HI,  and  which  incorporate  both  the  geometry  and  propagation  features  of  source, 
target,  and  the  medium  explicitly. 

Although  the  present  approach  includes  nongaussian  noise  components,  it  is  for  the 
moment  analytically  limited  to  iso-velocity  propagation  conditions,  where  Vc  =  0.  The  quantitative 
extension  to  the  more  realistic  geometries  implied  by  Vc  *  0,  namely  first-to-third  convergence 
zone  (CVZ)  applications,  is  one  of  the  suggested  next  steps  in  extending  the  analysis.  However, 
the  assumption  Vc  =  0  here  in  no  way  destroys  the  canonical  nature  of  the  results  or  the  generality 
of  the  methods  employed. 

*  127  East  91st  Street;  New  York,  NY  10128. 
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The  organization  of  this  paper  follows  closely  the  topical  outline  of  I-IV  above,  with  a  short 
summary  of  the  principal  results  and  needed  next  steps  in  the  analysis  and  suggested  experimental 
efforts. 


SECTION  2.  THE  CANONICAL  CHANNEL  AS  A  LINEAR  FILTER 
The  canonical  channel  is  shown  schematically  in  Figure  2.1.  It  is  postulated  to  be  linear, 
but  is  otherwise  unrestricted.  This  structure  suggests  that  it  can  be  represented  equivalently  by  a 
spatio-temporal  linear  filter  that  is  generally  time-  and  space-variable  with  a  variety  of  injected  noise 
sources,  in  the  manner  of  Figure  22  below,  where  hs(=  hR  +  hD)  is  the  weighting  or  green's 
function  of  the  medium,  excluding  the  coupling  operations  hT  Rc .  Here  hR  and  hD  are  respectively 

the  random  and  deterministic  components  of  hs  =  hs(R',R|x,t)  . 

The  quantities  hT,Rc(~*^T,R)  in  Figure  2.1  represent  the  coupling  to  the  medium 
(T  =  source;  Rc  =  receiver),  with  the  beam  patterns  associated  with  the  filters  (h-p^R)  of  these 

A 

coupling  apertures.  Here  also,  Sin  =  S|n(t)  is  the  original  (encoded)  input  signal,  and  GT  is  the 
aperture  operator,  such  that  (Sec.  II  [2]  the  input  signal  source  function  Gt,  becomes 

GT(t,y  =  GTS„  =  JT  hT(t-  tiysj?  ,  (2.1) 

where  in  general  hj  is  treated  as  time-variable,  because  of  possible  motion  of  the  source. 


1  a’ 


lat 


Inhomogeneous 

Medium 


b’ 


lAR 


(hRc,^R)dVR 

A 

=  R 


— Generalized  Channel-  — 


Figure  2.1  Operational  schematic  of  an  inhomogeneous  linear  medium,  with  source  and 
receiver  coupling.  Here  Q  is  the  inhomogeneity  operator  associated  with  this  general, 
inhomogeneous  linear  medium. 

A 

Similarly,  X(t)  is  the  output  of  the  receiving  aperture,  expressed  by  the  operator  R,  e.g., 

X(t)  =  Ra(R,t)  =  JVidVB(R)£  hRcl(t  —  x, tIR  e  VR)a(R,x)dx  (2.2) 
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where  hRC  (*hR  above)  is  possibly  time-variable  also.  Here  a(R,t)  is  the  field  propagated  in  the 
medium,  due  to  ambient  and  injected  signal  sources  (Sin).  Note  that  both  the  transmitting  and 
receiving  apertures,  or  arrays,  are  generally  frequency  selective:  only  for  suitably  narrowband 
signals  vis-k-vis  these  apertures  can  the  latter  be  treated  as  essentially  frequency  insensitive,  an 
approximation  used  throughout  this  study. 


2.1  Channel  Operators  and  Filter  Equivalents: 

The  various  channel  operators  of  Figure  2.1  are  represented  by 

Xui  =  R  =  —  GT  —  GA,  (2.3a) 

where  GA  is  associated  with  ambient  noise  in  the  medium  (Sec.  2.1  of  [1]),  and  that  of  the  medium 
itself  is  given  by  the  operator  relation  (Sec.  2,  [1]). 

Tm  ’  =  (L(0)  -  Q)”\  (2.3b) 

where  the  general  field  a(R,t)  is  formally  found  from 

(L(0)  - Q)a  =  -Gt  - Ga,  with  (L<0)  -  Q)g(Q)  =  -5tt,8RR„  (2.4) 

in  which  g(Q)  is  the  green's  function  of  the  medium,  including  boundaries  and  GA  is  the  local 
(space-time)  source  functions  of  the  ambient  field.  The  L(0>  -  Q  represents  the  propagation 
operators.  For  example,  for  an  inhomogeneous  medium  obeying  a  Helmholtz  wave  equation,  one 

A  A  A 

can  write  (in  rectangular  coordinates:  R  =  ixx  +  iyy  +  izz ): 


\ 


1  U 

V2 — 2-(l+e[R,t])— j  a  =  -(GT +  GA),  or  =  0,  outside  source  domains,  (2.5a) 
c„  ot  j 


with 


L<0)  =  V2  - 


i  a2 


1 


c?3t2  ;Q‘c2£(R,t)at2 


(2.5b) 


for  a  particular  representation  of  the  medium.  The  homogeneous  part  is  described  by  L(0) ,  while 

a 

the  inhomogeneous  effects  are  represented  by  Q .  Of  course,  the  medium  from  our  viewpoint  is 


random,  so  that  properly  Q  is  a  random  operator  (with  e  is  (2.5b)  a  random  field).  Equations 
(2.4),  (2.5)  embody  an  ensemble  of  representations,  which  is  called  a  Langevin  Equation,  whose 
solutions  are  statistics  of  the  random  field  a.  [See  [2],  Part  I;  also,  Chapter  10  of  [3].]  It  is  these 
statistics,  particularly  the  lower  order  moments,  which  we  shall  need  to  describe  the  effects  of  the 
generalized  channel  on  our  transmitted  and  received  signals  here,  contained  in  X(t),  (2.2). 
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Figure  22  Generalized  linear  inhomogeneous  channel,  showing  ambient  noise 
mechanisms,  scattering  operations,  and  coupling,  for  input  signal  Sjn  and  (temporal) 
receiver  processing  for  detection  and/or  estimation.  The  random  and  deterministic 
components  of  the  medium  .are  indicated. 


The  space-time  (now  random)  filter  equivalent  to  the  random  inhomogeneous  medium  is 
schematically  presented  in  Figure  22,  which  is  the  "filter”  form  of  Figure  2.1  above.  The  green's 
function  or  "impulse  response"  of  this  inhomogeneous  medium  is  indicated  by 

hs(R',RiT,t)  =  g<Q)(R,tlR',$)  =  g,05  +go 1 ,  (2.6) 

with  g(RQ),  g|)Q)  respectively  the  random  and  deterministic  parts  of  g(Q)  A  major  task  of  all 

investigations  of  such  media  is  to  determine  g(Q),  or  reasonable  approximations,  and  the  associated 
statistics.  The  results  of  this  must  then  be  applied  to  the  various  noise  (reverberation)  and  signal 
models,  from  which  in  turn  the  desired  detection  algorithms  and  system  performance  are  to  be 
obtained.  Unlike  the  purely  ambient  noise  cases,  the  dominant  noise  effects  are  signal-dependent 
here,  which  puts  a  premium  on  the  choice  of  transmitted  signal,  Sin. 

Before  we  proceed  to  results  for  the  dominant  scatter-channel  (Section  3,  et  seq.),  let  us 
observe  that  the  (linear)  space- time  filters  (~  g(Q))  presented  formally  so  far  can  be  replaced  by 
purely  (linear)  temporal  filter  equivalents,  under  often  reasonable  conditions.  This  can  be  useful 
when  we  seek  the  medium's  response,  e.g.,  green’s  function,  empirically,  and  attempt  to  model  it 
analytically  with  simple  functional  forms,  which  may  preserve  the  main  effects  of  the  medium  and 
still  yield  manageable  results  for  detector  design  and  performance.  In  [1],  Section  2.1,  it  is 
established  that  the  necessary  and  sufficient  condition's  are: 
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A .  Ambient  Sources:  The  necessary  and  sufficient  (n.+s.)  conditions  that  that  portion  of  the 

canonical  channel  (medium  +  coupling,  cf.  Figure  2.1)  which  contains 
ambient  sources,  local  and  distributed,  be  representable  solely  as  a 
(realizable)  temporal  filter,  is  that  these  ambient  sources  be  temporal  only 
and  independent  of  location,  in  the  source  domain  Aa>  e.g., 

GA(R,t)  =  GoA(t),cf.  Sec.  2B,[1] . 

B .  Scatter  Channels :  The  (n.+  s.)  condition  for  the  canonical  scatter,  i.e,  signal-dependent  noise 

channels  is,  like  the  case  of  ambient  sources,  that  the  original,  injected 
signals  that  excite  the  scattering  elements  be  independent  of  location,  cf. 
Section  2C,[1]. 

Specifically,  for  the  ambient  and  scatter  components  of  the  canonical  channel  shown  in 
Figures  2.1  and  2.2,  which  obey  A,  B  above,  the  equivalences  between  the  purely  temporal  and 
space-time  filter  weighting  functions  are: 

A.  Ambient  Sources  : 

(Eq.  (2.24), [1]) 

hA(t — x, t)  =  f dR J dR' J"  hR„(t-t'.t|R)g<Q>(R,x'|R',T)dt'.  (2.7) 

*  "a  •  aa 

B.  Scatter  Channels: 

(Eq.  (2.31a),[l]  \ 

hs(t  - x, t)  =  Jj dRdR/Jj  hRc(t -  x',  t|R|(Q)(R,x/|R,,t')hT(t/ -  x,x'|R')dx'dt' .  (2.8) 

Then,  for  the  output  of  the  receiving  aperture  or  array  (T^),  Figures  2.1  and  2.2,  where  now  the 
canonical  channel  is  described  solely  in  terms  of  purely  temporal  filters  (hA,  hs,  (2.7),  (2.8)),  we 
see  that  X(t)  is  given  by  the  (temporal)  convolutions. 

X(t)=  hs©  S^  +  hA©  GoA 

=  hs©  S^  -H  hd  o  Sk  +  NA(t)  (2.9) 

for  the  random  and  deterministic  portions  of  the  scatter  channel.  Note  that  NA(t)  can  have  a 
deterministic  part ,  e.g.,  NA  (t)  *  0  in  this  model.  Also,  the  ambient  noise  can,  in  pan,  be 
scattered  by  medium  inhomogeneities,  cf.  (2.7).  Physically,  this  corresponds  to  a  coherent  or 
specular  component  of  the  ambient  sources,  generated,  say,  by  reflections  off  the  ocean  bottom 
and/or  surface.  The  purely  temporal  equivalent  of  Figure  2.2  is  now 
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Figure  2.3  Purely  temporal  channel  equivalent  of  the  space-time  generalized  channel 
of  Figures  2.1  and  22. 

Special  cases  of  (2.9)  are  familiar  from  telecommunications  applications.  We  have 

non  -  gauss  stastics: 

gauss  statistics:  hs  or  g<Q)  (hs  or  g(Q) ) 

I.  hd=0;=  "Rayleigh  fading"  -  troposcatter,  Nakagami  models; 

ionospheric  chanlels:  unresolveable  same  physical  mechs. 

multipath  (scatter);  fading. 

II.  hd  *  0;=  "Ricean  fading"  channel;  VHF,  UHF,  Nakagami  models; 

HF  groundwaves,  etc;  hd(x,t)  =  hd(x).  DM  poisson  models. 

III.  hr  =  0;=  time  invariant  channel,  hs  =  8(x-t),  same  (2.10) 

plus  NA(t)  =  AWGN(=  add.  white  NA  =  "AWnonGN." 

gauss  noise);  hs(x),  purely  frequency  same 

selective. 


A  brief  look  at  the  implications  of  these  results  follows. 
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2.2 


Implications: 


A.  The  General  Space-Time  Formulation:  X  =  Ra,  Eqs.  (2.2) -(2.4), 

Figures  2.1  and  2.2: 

The  advantages  of  the  general  approach  are: 

1) .  Physical  completeness;  predictive,  from  propagation  models; 

2) .  Arrays  -  coupling  to  medium,  array  design,  etc; 

3) .  Spatial  as  well  as  temporal  effects;  explicit  inclusion  of  Vc  *  0,  scattering  etc.;  and 

4) .  Correct  propagation  of  ambient  sources  in  scattering  medium. 

The  major  disadvantage  is  the  difficulty  of  obtaining  g(Q)  (or  M(Q),  (3.2))  in  the  frequently 
encountered  heavy  scattering  situations  (e.g.,  multiple  scattering),  and  the  quantitative  treatment  of 
cases  Vc  ^  0  involving  a  number  of  convergence  zones.  These  limitations,  however,  are  already 
overcome  in  many  cases.  See  [1],  [6],  [8],  in  [1]  here. 

Further  progress  may  be  expected,  computationally  ([10],  [11],  in  [1]  here),  as  well  as 
analytically. 

B.  The  Purely  Temporal  Formulation:  X  =  hs©  Sln  +  NA  ,  (2.9): 

The  principal  advantages  of  this  more  restricted  channel  model  are: 

1) .  Analytical  and  (comparative)  statistical  simplicity; 

2) .  Ad  hoc  choices  of  hs  and  the  associated  statistics  justified  by  "local"  experiments,  i.e., 

measurements  made  at  the  receiver,  for  specific  channels; 

3) .  Direct  ties-in  to  many  of  the  channels  and  "target"  modeling  results  of  the  1960's  and 

1970's  for  radar,  sonar,  telecommunications,  etc.,  Chapter  2,  [4],  Chapter  13,  [5]. 
For  example,  a  simple  often-used  model  is  the  wide  sense  stationary  uncorrelated 
scatter  (WSSUS)  channel,  cf.  Sec.  3  ff. 

The  main  disadvantages  are: 

1) .  Tenuous  relations,  at  best,  to  the  underlying  physical  mechanisms:  no  insights  as  to 

what  is  happening; 

2) .  Dependence  on  particular  channel  measurements  -  non-predictive  for  others; 

3) .  No  explicit  spatial  structure:  does  not  permit  array  design  or,  more  broadly,  (local) 

matched  field  processing*  (the  spatial  equivalent  of  "matched  filtering"  in  time). 


$ 

This  is  not  the  same  as  "Matched  Field  Processing"  (MFP),  as  used  in  passive  source 
localization;  see,  for  example,  [18],  [19]  of  [1]. 
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3.  SCATTER  CHANNELS:  A  SUMMARY  OF  RESULTS  [1],  [2] 


The  received  wave  X(t),  after  a  signal  Sto(t)  has  been  injected  into  the  general  channel,  cf. 
Figure  2.1,  has  been  represented  by  Eq.  (2.9).  The  space-time  equivalent  can  be  expressed 
compactly  in  operator  form  by 

X(t)  =  M(Q)GT(-Sin)  +  RM(Q)(-Ga)  ,  (3.1) 

where  M(Q)  is  the  integral  green's  function  associated  with  (2.4),  viz, 

M(Q)(R,tlR',t')  =  -JJ dt'dRyQ)(R,tlR',t')(  )R,  t.  .  (3.2) 

This  operator  is  formally  described  by 

M(Q)  =  (l  -  f|)-1  ;  tj  =  M^Q  ,  (3.3) 

where  Mm  is  the  integral  green's  operator  of  the  associated  homogeneous,  i.e.,  infinite,  medium, 
in  which  the  inhomogeneity  operator,  Q,  vanishes,  i.e.,  g(Q)  — >  gM  from  (2.4)  with  Q  =  0.  Thus, 
is  defined  by  Eq.  (3.2),  with  now  g{Q)  — >  g^.  These  relations  hold  formally  for  random 
fields,  as  well,  cf.  remarks  following  Eq.  (2.5b),  where  g(Q),  f),  M(Q>  are  stochastic  functions  and 

operators;  is  deterministic.  The  full  field,  a(R,t),  including  both  the  "homogeneous"  or 
unscattered  part,  aH,  and  the  "inhomogeneous"  scattered  component,  ax  =  a  -  aH,  is  formally 
obtained  as  the  solution  to  Eq.  (2.4),  namely, 

a(R,t)  =  (l-fi)  1aH  =  ^  Tj(n)aH;  aH(R,t)  =  Moe(-Gx) .  (3.4) 

n=0 

Equation  (3.4)  is  the  Langevin  (ensemble)  equation  for  the  random  field  a  ,  and  as  explained 
above  (Section  2),  its  "solutions"  now  are  the  various  statistics  of  a  ,  viz.,  (a),  (a2),  (a,a2), 

Wj  (g),  etc.  These  are  the  quantities  needed  in  determining  effective  models  of  the  desired  (target) 
signal  and  the  accompanying  reverberation  and  ambient  noise,  in  order  to  approach  optimal,  or 
practical  near-optimal,  detection  and  classification,  as  noted  in  more  detail  in  Section  5  ff. 

3.1  SECOND  MOMENTS  AND  WSSUS  CHANNELS 

Unfortunately,  fully  analytic  solutions  for  g(Q),  M(Q) ,  and  therefore  a  and  its  statistics,  are 
not  generally  achievable.  We  must  accordingly  seek  compromise  and  approximate  solutions.  In 
the  rest  of  this  section,  we  summarize  some  of  the  principal  results.  We  are  concerned  principally 
here  with  narrowband  injected  signals.  For  the  details,  see  Section  3  of  [1]. 

A  statistic  of  major  importance  in  describing  the  effects  of  our  canonical  channel  on  the 
injected  signal,  S^t) ,  is  the  cospreading  function  Gs,  defined  here  by 


8 


(3.5) 


QsCtj.Xj;  Vi)Hs(x2,  v2)*)  ,  (Eq.  (3.3c),  [1]) , 

where  Hs(x,  v)  is  the  well-known  spreading  function  ([4],  [5]),  which  measures  how  much  the 
channel  delay  spreads  in  time  (~  x)  and  spreads  in  frequency  (~  v)  whatever  is  transmitted  through 
it 

When  multiple  scatter  can  be  ignored,  i.e.,  only  single-scatter  is  the  dominant  component 
of  the  reverberation,  then  Eq.  (3.4)  becomes  the  (Bom-)  approximation 

a(R,t)  =  aH  +M.QaH  »  with  a,(R,t)  =  M_QaH  (3.6) 

the  reverberation  contribution.  Then  it  can  be  shown  (Section  3.2,  [1],  and  references)  that  the 
cospreading  function  Eq.  (3.5)  takes  the  form 

Qs.(Xi,x2;  vl,v2)  =  as(x„v1)8(x2-x,)5(v2-v1) ,  cs>0,  (3.7) 

which  is  the  result  of  independent  scattering,  namely,  scattering  with  independent  emissions,  path 
delays,  locations,  and  Dopplers.  The  quantity  as  is  usually  called  the  scattering  Junction  and  with 
(3.7)  defines  completely  the  second-order  moment  statistics  of  the  channel  (cf.  refs,  in  [1]). 
Channels  for  which  (3.7)  holds  are  called  WSSUS  or  Wide  Sense  Stationary  Uncorrelated  Scatter 
channels  (cf.  p.  18;  [4]),  and  are  often  used  to  model  EM  telecommunication  and  underwater 
acoustic  channels  (for  example,  V,  [2];  Sec.  2.  [6]).  Derivation  of  the  scattering  function  is  one 
major  goal  of  scattering  analysis,  cf.  Sec.  3.2  following.  However,  WSSUS  channels  do  not 
usually  appear  to  be  good  models  over  long  ranges  in  the  ocean  [7],  and  may  not  be  for  the  shorter 
ranges  considered  here.  The  latter  needs  to  be  investigated,  particularly  when  multiple  scattering 
may  be  significant;  see  the  remarks  in  Secs.  3.3,  ff. 

An  important  result  in  the  WSSUS  cases  is  the  covariance  of  the  received  wave, 

Kx(x,,x2)  =  (XjX2)  -  (X1)<X2> ,  from  which  the  intensity  and  power  spectrum  can  be  directly 
obtained.  For  narrowband  signals,  this  WSSUS  covariance  and  the  associated  scattering  functions 


are  related  by  (cf.  Secs.  3.2-3.4,[l]) 

=  ExRe|eto,(Xl  Tj)JJors(t’v)(so(ti  ~x)ia 

^e2*v(T1-T1)dtdvJ>(3.8) 

=  ExRc{kl1)(x1,x2)exp(ico0[x2  -x,])} , 

(3.8a) 

where  s„  is  the  complex  envelope  of  the  initially  injected  narrowband  signal*  (about  the  "carrier" 
/0).  Here  x  ~  Toj  ,  the  path  delay  from  source  to  scatterer  to  receiver,  in  the  far- field,  and  Ex  is 
the  mean  intensity  of  the  array  (or  aperture)  output,  X(t) ,  related  to  cts  by  as  =  as  /  Ex  ,  where 


*  "Narrowband  "  is  defined  here  as  a  ratio  of  central  frequency  f0  to  bandwidth  A/ ,  which  obeys  f0/Af  >  4, 
approximately. 
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Ex  =  JJ °S(T, v)|s0(tj  -  T)J2dT  dv  .  (3.8b) 

The  average  (  )  is  an  ensemble  average  over  the  phase  or  epochs  of  the  injected  signal 
envelopes  (see  Eqs.  3.10  and  3.11  ff).  In  terms  of  the  ambiguity  function ,  defined  here  by 

X>v> s  -  D-CAaV  +  ,<o  =  2  rf ,  (3.9) 

where  X.(o.o)=f.k  (t\\  dt'  =  1  is  the  intensity  of  the  normalized  envelope  amplitude  s0  ,  we 
can  express  the  WSSUS  covariance  (3.8)  as 

K?(t1,tJ)  =  Exjss(-v,tJ-t,)X0(t1-t„v)Vi”iv,Mv  (3.10) 

Here 

Ss(f,t) «  F>",f,{as(t.v)}  =  JJ  os(t,  v)e-iwt+2*ivtdT  dv  (3.10a) 

is  the  double  Fourier  transform  (f )  of  as ,  as  indicated. 


3.2  FULL  CHANNEL  COVARIANCE:  NARROWBAND  SIGNALS 

One  of  our  main  tasks  in  channel  modeling  is  to  determine  not  only  the  WSSUS  covariance 
(Eqs.  (3.8)  and  (3.10)),  but  the  complete  covariance,  including  the  contributions  of  multiple 
scatter.  This  can  be  noticeably  advanced  with  the  help  of  a  suitable  counting  functional ,  which 
allows  us  to  specify  the  (formal)  interactions  of  individual  scattering  elements  with  one  another  and 
the  incident  radiation.  The  detailed  development  is  provided  in  Secs.  3.2  and  3.3  of  [1].  The 
complete  covariance  is  found  to  be  (cf.  (3.30), [1])  for  these  narrowband  injected  signals 


Kx(ti.t2)=2Re1 


i<°°(ti  t^JJ^s(^i.v'i)e27nVlft2  ^(Mti  -  x2)*inSj^dtldvl 

[(so(tl  -  *1)^80(12  -  T2)^)  dTj...dv2 


(3.11) 


v.  J 

The  first  term  embodies  the  WSSUS  or  single-scatter  component,  while  the  second 

oo 

represents  all  the  multiple-scatter  contributions  (k  >  2)  ,  viz.  >  where  is  now  the 

k-2 

cospreading  function,  cf.  (3.5),  for  the  latter  interactions. 

Using  our  earlier  results,  cf.  Eq.  (4.47)  of  [8],  for  narrowband  (n.b.)  signals  in  the  far- 
field  (f.f.),  we  have  for  a  typical  bistatic  received  waveform  from  the  j^1  point  scattering  element 
here 
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u 


.b=  ,^f"^-exp[i(0.  (1  -  e)(t  -  T„  a))]^T^„ 
(4tc)  RtRr 


n 

f.f. 


so(t-T0(X)y 


(3.12) 


l/o 


where  e  =  (~  vd  /  c0)  is  a  random  Doppler,  |i  =  1  -  e ;  yo(X)  =  "cross-section"  of  the  scatterer, 
and  AT,AR  are  the  (complex)  beam  patterns  (cf.  Figure  2.1 )  of  transmitter  and  receiver,  e.g., 

At  =  Ar(vT  -  voTl /„)  ,  with  vT  =  iT/„  /c0  ,  and  voT  is  a  steering  wave  number,  etc.  For  details 
of  the  geometry,  see  Figure  1  of  [2].  Here  X  is  a  distance  (=  |R|  /  c0)  measured  in  units  of  time. 

Thus,  T0(X)  =  (|Rt|  +  |Rr|)  /  c0 .  For  use  in  the  time-frequency  forms  above,  Eq.  (3.5)  et  seq.,  we 
can  replace  T„  by  t ,  co0  by  2tcv  ,  etc.  (cf.  comments  following  Eq.  (3.28),  [1]).  From  Eq. 

(3.12),  then,  it  is  readily  shown  that  for  the  WSSUS  component,  the  scattering  function  cs  now 
becomes  explicitly 


A?  1- 


2tcv 


as(t,v)  = 


CO 


0  / 


ArAj)  G(t) 

/T,V 


>0 , k=l , 


(3.13) 


(47t)4RT(x)2RR(t)2 

with  G(x)  the  contribution  of  the  Jacobian  |J|  from  rectangular  to  spherical  coordinates  [8],  and  the 
conversion  of  X  =  X(x),  dX  =  dx  /  f '(x) ,  etc.,  all  of  which  depend,  of  course,  on  the  geometry  of 
source,  receiver,  and  the  illuminated  scatterers  on  the  boundaries  and  in  the  volume  [6],[8]. 

Various  useful  generalization  and  special  results,  which  are  beyond  the  scope  of  the  present 
summary,  as  well  as  extensions  of  the  methods  employed  to  derive  Eqs.  (3.8)  and  (3.1 1)  above, 
are  presented  in  detail  in  Section  3.3  and  Section  4  of  [1]. 


3.3  SUMMARY  OF  PRINCIPAL  RESULTS  FOR  CHANNEL 
MODELING:  SECS.  2,  3,  AND  SECS.  2-4  OF  [1] 

To  operate  effectively  in  a  physical  channel  environment  (=  medium  plus  coupling),  we 
must  estimate  the  pertinent  features  of  the  channel.  Thus,  our  ultimate  signal  processing  is 

necessarily  adaptive.  Accordingly,  probing  or  "remote  sensing"  of  the  channel’s  response  (hs)  is 
suggested.  As  long  as  the  channel  is  stable  during  intervals  of  use,  we  may  expect  good  results. 
Over  time,  we  require  updating  of  channel  parameters,  which  is  usually  done  on  a  continuous 
averaging  basis.  Without  adaptivity,  we  cannot  make  full  use  of  our  processing  capabilities.  We 
may  then  expect  serious  degradation  of  performance  vis-a-vis  the  adaptive  cases,  which  in  the 
broad  sense  are  "matched”  to  actual  channel  characteristics.  In  the  above  when  purely  temporal 
specification  of  the  channel  is  employed,  through  the  channel’s  response  function  hs(t,x) ,  we 
have,  in  effect,  preformed  beams:  the  array  or  aperture  has  a  postulated  structure,  so  that  only 
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"matched  filtering,"  i.e.,  temporal  matching  for  optimization,  is  possible.  On  the  other  hand,  if  the 
coupling  structure  is  left  adjustable  to  the  incoming  signal  and  noise  fields,  then  spatial  matched 
filtering,  or  (local)  "matched  field"  processing*  is  possible  as  well.  For  this,  of  course,  we  need 
pertinent  statistics  of  the  noise  (and  signal)  fields,  such  as  their  covariance  functions. 

It  is  with  this  background  that  we  summarize  our  principal  results  above  on  channel 
modeling,  which  may  serve  as  a  guide  to  both  modeling  and  measurement,  for  ultimate  application 
in  the  required  signal  processing.  We  note  the  following  : 

I .  For  there  to  be  a  purely  temporal  representation  of  the  channel  equivalent  to  the 
fundamental  space-time  representation,  the  driving  sources  (ambient  and  active)  must 
be  independent  of  location.  This  is  both  necessary  and  sufficient,  cf.  Sec.  2.1, 

A,  B. 

II.  The  temporal  channel  equivalent  (hs(t,x))  is  "global"- contains  all  spatial  effects- 

while  the  green's  function,  gCQ)  ,  Eqs.(2.4),(2.6)  of  the  medium  itself  is  "local," 
i.e.,  at  a  point  P(R,t) ,  although  its  value  at  a  point  includes  the  effects,  at  that  point, 
of  all  scatter  interactions  elsewhere  as  well. 

III.  Scattering  media  (including  boundaries)  are  nonreciprocal :  source  and  receiver  are 
not  interchangeable,  and  no  generalized  Huyghens  principal  (GHP)  exists  for  the 
solution  of  the  field  propagation  equations  in  such  media.  Alternative  (operational) 
methods  must  be  employed  [9],  cf.  Eq.  (3.6)  of  [1]. 

IV.  Single-scatter  channels  are  always  WSSUS  (in  the  steady  state).  They  accordingly 

have  a  cospreading  function  ,  which  obeys  Eq.  (3.7)  and  defines  a  scattering 

function  as(x,v) ;  (see  Refs.  [10]  and  [11]  for  such  a  surface  scatter  mechanism 
generated  by  wind). 

V.  WSSUS  channels  may  not  be  satisfactory  models  for  long-range  (to  CVZ  and 
beyond)  underwater  acoustic  propagation ,  because  of  coupling  (multiple-scatter 
effects,  which  are  path-dependent. 

VI.  Accordingly,  the  WSSUS  character  of  the  channel  should  be  empirically  tested~as 
part  of  the  adaptive,  remote  sensing  for  channel  matching;  (see  Refs.  [10]  and  [  1 1] 
for  surface  scatter). 

VII.  The  scattering  function  as(i,v)  can  be  calculated,  as  well  as  measured;  (see  Sec.  3.2 
above). 


* 


See  footnote,  p.  7,  for  comments. 
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VIII.  The  general  second-order  (output)  covariance  Kx(tj,t2)  is  explicitly  obtained 
(Sec.  3.2)  when  there  is  a  non-WSSUS,  as  well  as  WSSUS,  component. 

IX.  Generalizations  of  VIII  are  derived  in  Sec.  3.3  of  [1]. 

X .  Extensions  and  equivalences  with  counting  functional  results  are  derived  in 
Sec.  4  [1]. 

Finally,  the  above  suggests  some  further  topics  for  study  along  these  lines: 

1 .  Extension  of  WSSUS  concepts  to  include  spatial  variations:  (local)  "matched  field" 
concepts. 

2.  Extension  of  the  global  channel  representation  (~  hs)  to  specific  evaluation  of  media: 
media  "probes." 

3.  Further  direct  calculation  of  specific  classes  of  scattering  function  as ,  and 
cospreading  functions  (based  on  [2],  [6],  and  [8]),  as  well  as  on  the  approach 
outlined  in  detail  in  Sec.  3.2,  [1]. 


4.  ELEMENTARY  TARGET  MODELS:  TARGET  STRUCTURE  FACTORS 

(SECS.  6,7  OF  [1]) 

The  aim  is  to  construct  simple  but  "realistic  target  models,  which  exhibit  reflection  and 
,  scattering.  These  models  are  "quasi-phenomenological"  in  some  details,  directly  physical  in 

others.  The  basic  idea  is  to  preserve  the  relevant  target  features,  e.g.,  "cross-section,"  spatial 
distribution,  motion,  which  are  dominant,  as  is  the  central  role  played  by  the  associated  geometry. 
Much  of  the  so-called  "fire-structure"  is  smoothed  out  by  the  inevitable  uncertainty  in  parameter 
values,  through  statistical  averaging.  This  is  not  a  new  problem,  nor  are  our  results  here 
particularly  innovative:  the  comparatively  naive  model  employed  here  consists  of  a  linear 
distribution  of  reflecting  "sources,"  with  adjustable  parameters  of  scale  and  function. 

!  Our  particular  goal  is  to  obtain  analytic  expressions  for  the  received  narrowband  (i.e., 

complex)  envelope  (s0  =)sQ(t)L  ,  Eq.  (3.8)  et  seq.,  of  a  linear  target,  arbitrarily  oriented  with 
respect  to  the  transmitter  and  receiver,  in  the  manner  of  Figures  4.1  and  42. 
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(This  is  Figure  1  of  [6],  with  the  surface  area  replaced  by  a  line,  L.) 

K 

I 


I 


Figure  42  Target  orientation  vis-k-vis  Op,-*  Op,  Or  . 


The  result  (cf.  Eq.  (6. 19)  of  [1])  is  the  continuous  linear  model,  for  which 


(4.1) 

is  the  complex  signal  envelope  of  the  target  as  seen  at  the  receiver.  Here  *4.^  =  ARAT  ,  while 


>.  ->  s„(t-T„;e)L  =exp{i0),,(To-e)}^G<i>(X)|L^ra(<)e^ik'‘c‘"-' 


GW(/“)HfSrS?’Eq- 


(4.1a) 


*  A  A 

and  £  =  it£=  i Xlzlx  +  iytxl2  +  A.2.3  the  direction  cosines  of  l  vis-^-vis  (x,y,z)  of  0T  . 

The  quantity  t  in  Eq.  (4.1)  represents  different  possible  Doppler  effects,  specified  by  Eq.  (6.10)  of 
[1],  Other  quantities  in  Eq.  (4.1)  are  given  by  Eq.  (4.5)  ff: 


(4.2a) 
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vD  =  doppler  velocities:  for  no  platform  motion, 

Vjjt  -4  (t  -  Ror/cJv,,  ;  but  see  Eq.  (67.10b),[l] ;  (4.2b) 

To  =  (Rot  '  ^qr)/Co 

=  (mean)  wavefront  path  time:  Ot  -4  0L  — »  Or  ;  (4.2c) 

» 

c0  =  speed  of  sound  in  (isovelocity  medium  (Vc  =  0)) ;  (4.2d) 

s  =  27cif ;  s0  =  27tif0 ; 

fc  =  central  frequency  of  injected  (n.b.)  signal;  (4.2e) 

a  =  absorption  coefficient;  (4.2f) 

vot»vor  =  (wave  no.)  beam  steering ;  (4.2g) 

=  ( preformed ')  beam  patterns  (which  are  complex) ;  (4.2h) 

(See  also  Appendix  I,  [6],  for  definitions  and  relations), 

e.g.,  t'  =  t- AtdS,  cf.  (3.14),  [6].  (4.2i) 

The  quantity  h0  is  the  (real)  narrowband  portion  of  the  (real)  target's  scatter  kernel 

qL(x,t)  =  qL(t-T)  =  2h0(t-x)LRe{exp[ico0(t-T)-iy0(t-T)]}  ,  (4.3) 

cfEqs.  (6.9),  (6.18)  of  [1]. 

An  alternative  to  the  continuous  aperture  or  scatter  kernel  qL  (4.3)  is  the  discrete  model , 
where  now  the  aperture  takes  the  form 

qL  (x,  t)  =  qL (t  -  x)  =  X a kei$k  5(^  ~  i*kA 4)8 (t -  X  -  T0 ) ,  (4.4) 

k>l 


A  A 

with  i,  =  unit  vector  along  t  as  before,  cf.  (4.1)  et.seq.,  and  Figure  4.2,  e.g.,  i  =  i/  above.  The 
discrete  equivalent  of  (4.1)  is  then  found  to  be  (cf.  sec.  7,  [1]) 


S„(t-T„  ;  e)t  *  ^G'e“,‘‘X^A„(kAey2,a- 

k-1 

•exp[ik0(2a0  i/)vdRoR/c0  so(t-T0-e)in  , 


(4.5) 


since  =  ivd  .  In  the  far-field,  *4RT(kA^)  -4  ,  a  constant  (dependent  on  the  "canier''  or 

central  frequency,  f0  and  target  length  L).  Here  the  weights  ak  =  (ak  +  Aak)  may  be  regarded  as 

mixed  deterministic  and  random,  at  positions  t  =  kA^  ,  k  =  1 . K  along  L  .  For  example,  the 

Aak  may  be  regarded  as  independent  gaussian  quantities,  or  the  ak  ,  each  k  ,  may  be  regarded  as 
obeying  Poisson  statistics 
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with 


Wi(ak)  = 


e  k  ,  A£vk  =ak  ,  etc., 


(4.6a) 


m! 


4>k:wi($k) »  .  (4.6b) 

to  be  selected,  appropriate  to  the  model  chosen. 

The  quantity  F0(t)  is  called  the  Target  Structure  Factor  (TSF)  and  is  defined  by  the  ratio 

(4-7) 


for  these  coherent  reception  modes,  with  sc.L  given  specifically  by  Eqs.  (4.1)  or  (4.5)  here.  As 
we  shall  see  presently  (Secs.  5.2.A,B,  Fc  plays  a  central  role  in  effective  detection  in  a 
reverberation  dominated  scenario. 

Various  extensions  of  the  models  (4.1),  (4.5)  need  to  be  considered.  Not  necessarily  in 
order  of  importance,  among  these  are: 


1.  Moving  target  and  platform  scenarios  ; 

2.  Deterministic  and  random  dopplers ; 

3.  Explicit  statistical  models  of  scattering  strengths : 
h0 ,  (4.1),  ak  ,(4.4); 

4.  Intensities  and  covariances  of  these  generic  target  models 
for  incoherent  detection ; 

5  Continuous  models  (~  h0.L  ,  (4.1)) ,  (4.3)  continuous  in  i{e  L)  ; 

6.  More  detailed  physics  of  qL  ,  (4.3);  "highlights,”  wakes,  resonances 
etc.  ;  and 

7.  Kirchoff  boundary  conditions  on  l . 

For  a  more  complete  account  of  these  models  and  derivations,  see  Secs.  6  and  7  of  [1].) 


5.  TARGET  DETECTION:  THRESHOLD  ALGORITHMS  AND 
PERFORMANCE  MEASURES 

Our  aim  in  this  section  is  to  provide  a  concise  summary  of  analytical  results  for  optimal  and 
near-optimal  detection  of  target  signals  in  noise  and  interference.  Specifically,  we  are  concerned 
with  the  combination  of  ambient  (gauss  and  nongaussian)  background  noise  and  signal-dependent 
noise,  or  reverberation,  when  the  latter  is  strong  vis-^-vis  the  desired  target  signal.  The  principal 
noise  or  interference  here  is  canonically  nongaussian,  in  the  sense  of  the  author's  familiar  Class  A 
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5.1  THRESHOLD  DETECTION  IN  SIGNAL-DEPENDENT  NOISE 

We  begin  with  a  short  summary  of  the  detection  situation  in  signal-dependent  noise, 
represented  by  the  usual  binary  "on-off"  hypothesis  models  here.  These  are 


I.  General  On-Off  Detection:  Signal-Dependent  Noise: 

This  is  expressed  by: 

H0:  (no  signal):  N(S)  +  NA  +  NR  vs.  H,:  (signal  +  noise):  N(S)  +  NA  +  NR  +  S. 

t  U? 

reverb,  ambient  receiver  (5.1) 

noise 


H.  Strong  Scatter:  Approximate  Decision  Situation: 


H0:N(S)  +  Na+g  vs.  Hi:N(S)  +  Na+g  +  S  , 


(5.1a) 


where  NA+G  =  NA  +  N0  is  a  mixture  of  gauss  (N0)  and  nongaussian  noise  in  the  usual  way,  with  a 
usually  dominant  nongaussian  component.  The  reverberation  [N(S)]  is  either  gaussian  or 
nongaussian,  depending  on  geometry  and  ("carrier")  frequency  of  these  n.b.  signals  [15]. 

•  In  general,  the  intensity  of  the  accompanying  noise  can  be  represented  here  by 

=  'P(S)  +  ^A+G  =  N(S)2  +  N2a+o  .  (5.2) 

where 


¥(S)  =  N(S)2  =  (Fs  +Fy  +  Fb)AJ,  F  >  0  .  (5.2a) 

Here  Fs,  Fv,  FB  are  respectively  the  "structure  factors"  of  the  acoustic  reverberation  off  the  ocean 
wave  surface,  the  volume,  and  the  bottom.  At  least  two  of  these  terms  (including  Fv)  will  be 
present,  again  depending  on  the  geometry  of  source  (cf.  Figure  5.1,  the  possible  target,  and  the 
dispersive  character  of  the  medium.  For  many  of  our  applications  the  surface  component 

dominates,  so  that  'F(S)  =  A„  Fs  +  ^A+G  .The  detailed  nature  of  these  structure  factors  is 
determined,  likewise,  by  geometry,  source  frequency,  and  waveform,  and  the  appropriate  scatter 
model  which  governs  these  reverberation  terms;  see  Section  5.3  ff.  and  references. 

Figure  5.1  shows  some  typical  propagation  geometries,  where  Vc  *  0;  (here  c  =  c„  +  Vc(R), 
is  the  speed  of  wavefront  propagation).  Depending  on  target  location,  the  various  components  of 
reverberation  may  be  expected  to  occur.  We  remark  that  bottom  reverberation  is  almost  always 
nongaussian,  while  surface  and  volume  reverberation  can  be  gaussian,  if  the  effective  number  of 
ensonified  scatterers  is  sufficiently  large  and  no  few  scatterers  predominate  over  the  others.  At 
high  frequencies  and  small  angles  (primarily  for  surface  and  bottom  scatter),  however,  we  may 
frequently  expect  a  small  number  of  large  scattering  elements  ("facets")  to  dominate  in  the  beam,  so 
that  an  often  strongly  nongaussian  reverberation  results. 
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Figure  5.1  Sketch  of  some  typical  ("high-frequency"  or  n.  b.)  propagation  of 

beamed  signals,  when  Vc  *  0,  showing  potential  mechanisms  of  reverberation: 
T,  R  =  transmitter,  receiver,  Tg  =  target. 


A.  On-Off  Binary  Detection  Algorithm;  Coherent  Threshold  Detection  -  Canonical 
Forms: 

From  earlier  work  [2]  the  canonical  form  of  the  optimum  "on-off  threshold  detection 
algorithm  is  well-known  to  be  [2],  [14]. 


M,N 


log  Aj(xl0)  =  gj  (x)^  =  b£2*  +  log  |i- 

5Xxra.„)(em.n) , 

m,n*l 

(5.3) 

1  *  2 

Bi-coh  bias  term’  (space)  m  =  1,...,M;  (time)  n  =  1,...,N;  |X  =  p  /  q; 

(5.4a) 

6J-coh3varH.gSl=2L(*)J 

{j±<^>)2<sr>2}. 

v  tn  ,n  j 

(J  =  MN) 

(5.4b) 

i 

i 

^  ®  ^ min-coh 

(5.4c) 

i 

i 

"processing  gain” 

"min.  detectable 
signal  (power) 

)H  =  £1  /(^(xfHoJdx  ; 

(5.4d) 
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(5.4c) 


(0.,.)  =  (aira,)(‘.Xs'm’('.)).  •  "i*  (4m|2)  =  1  ■ 

(cf.  remarks  ff.  Eq.  (60d),  [2],). 

with 

-  a?1  ■  <*.  (5-12b)ff  (5.40 

in  which,  as  before  [2],  [13],  m  denotes  the  m^1  sensor  or  "space-sample,"  while  n  represents  tn 
(=nAt),  the  n^1  "time-sample,"  in  the  space-time  observation  interval.  Here  we  have 

x  =  X/V2T7  =  {xm  „}  =  {xj  ,  with  'F  =  Eq.  (5.2);  j  =  m,n,  etc.,  (5.5) 

where  x  is  the  normalized,  sampled,  input  field  data.  Also,  p  (=p/q)  is  the  ratio  of  a  priori 
probabilities  that  the  input  does  or  does  not  contain  a  desired  signal  (target).  Finally,  w,(x|H0)  is 

the  (first-order)  pdf  of  x,  the  typical  noise  data  sample,  which  usually  consists  of  an  additive 
mixture  of  nongauss  (reverberation  and  ambient)  noise  and  a  normal  component. 

Figure  5 2  shows  the  flow  diagram  of  the  optimum  threshold  algorithm  g* ,  here  for 

coherent  detection  in  generalized,  e.g.,  signal-dependent,  nongaussian  noise.  As  usual,  coherent 
detection  requires  that  the  signal  epoch  e  =  Eq  ,  a  known  quantity  at  the  receiver,  adjusted  to 

maximize  the  carrier  waveform,  [13],  [14],  [4].  The  various  components  of  g*  represent  the 
following  signal  processing  operations: 


jsK:H,} 


Figure  5 2  Flow  diagram  of  the  canonically  optimum,  coherent  threshold  detection 
algorithms,  with  adaptive  beam  forming  (m),  generalized,  signal-dependent  noise. 
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Eq.  (5.6): 

A 

(1)  R  =  array  operator:  samples  the  input  acoustic  field  a(R,t),  at  points 

R  — >  { rm } ,  m=l,...,M;  see  Part  II,  Section  3.1  of  [1]. 

(2)  (  )n  =  time  sampling:  the  result  is  a  (normalized)  space-time  sampled  waveform 

sample,  which  then  undergoes  the  signal  processing  described  below: 

(3)  -/(ZMNL)  =  "matching-to-the-noise,"  via  /(x),  Eq.  (5.4d),  to  yield  yj=m  cf. 

©  in  Figure  5.2. 

(4)  ®(6mfn)  =  cross-correlation  with  appropriate  model  of  the  received  signal,  i.e.,  one 

that  includes  such  effects  as  possible  "Doppler  smear,"  "fading"  or 
unresolvable  multipath,  and  beam  steering,  etc.,  cf.  Sec.  5.2A,  embodied 
in  the  Target  Structure  Factor  (TSF),  cf.  (5.11),  (5.7). 

(5)  Em  =  "beam  forming,"  when  combined  with  appropriate  path  delays  (~rm/c0)  in 

A  A 

R  (2.2)  to  produce  a  temporal  wave,  yn.  In  conjunction  with  R  and 
proper  delays  (in  <0m,n>),  cf.  ©  in  Figure  52,  this  represents  adaptive 
beam-forming. 

(6)  Zn  =  in  conjunction  with  (4)  constitutes  "matched-filtering,"  or  "matching-to-the 

signal,"  cf.  ©  in  Figure  5.2. 

(7)  ©Bias  =  addition  of  a  suitable  bias,  B* ,  to  ensure  both  optimality  and  consistency 

of  the  test  vs.  H0 . 

(8)  g*  =  a  number,  produced  by  the  operations  (1)  -  (7),  which  is  then  compared  to 

a  threshold  K  which  is  preselected  to  yield  a  predetermined  false  alarm 
probability,  ap.,  cf.  Sec.  5.2  ff.  The  resulting  decision  is  H0:  g*  <  K  or 
Ht:  g  >  K ,  etc. 

We  remark  that  the  system  produces  decisions  continuously  as  the  input  data  stream  is  continuously 
updated,  for  the  given  (moving)  period  (to,t0  +  T);  or  one  can  operate  it  in  a  series  of  decisions  at 
intervals  T,  etc. 

More  important  is  the  critical  feature  of  all  these  optimum  (and  similarly  related 
suboptimum)  systems,  namely,  that  they  are  adaptive:  they  require  the  measurement  and  updating 
of  the  relevant  signal  and  noise  parameters.  These  are,  for  example,  noise  level,  received  signal 
level  (~Ao),  the  (nongaussian)  noise  parameters,  PA  B,  associated  with  the  governing  class  of 

noise,  expressed  in  terms  of  w^xIHq)  here,  cf.  [14],  whether  it  be  Class  A  or  Class  B  -  probably 

Class  A,  for  a  dominating  reverberation.  (For  details  see  [14];  also,  Part  II  of  [2]  and  refs.)  In 
essence,  the  effectiveness  of  these  coherent  (and  incoherent)  detection  algorithms  depends  on  how 
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well  we  can  estimate  the  relevant  parameters  of  the  "environment"  so  that  good  adaptation  is 
achieved.  Provided  the  time-scale  of  change  of  the  environmental  parameters  is  slow  (i.e., 
"isothermal  rates")  vis-a-vis  the  observation  period  needed  to  achieve  effective  levels  of  correct 
signal  detection  probabilities  —  which  is  fortunately  the  case,  most  of  the  time  —  our  threshold 
systems  here  can  practically  approach  optimality  quite  closely  [13]. 

Finally,  we  note  that  those  optimum  threshold  systems  no  longer  remain  optimal  for 
stronger  (to  strong)  signals.  However,  on  an  absolute  basis  they  are  quite  satisfactory,  as  long  as 
care  is  taken  to  avoid  possible  destruction  of  the  information  in  the  desired  signal,  which  can  occur 
because  of  the  saturating  nature  of  the  (ZMNL)  /(x).  In  general,  larger  detection  probabilities  result 
for  increased  signal  strength  under  these  conditions. 


5.2  Performance  Measures 

In  order  to  predict  the  effectiveness  of  our  detection  algorithms,  (5.3)  et  seq.  and 
Figure  52,  we  need  suitable  performance  measures.  These  are,  comprehensively  for  most 
purposes,  the  probability  of  correctly  detecting  the  presence  of  the  (target)  signal  (S).  For  these 
threshold  cases  (which  likewise  provide  algorithms  for  the  stronger  signals,  cf.  end  of  Sec.  5.1 
above),  it  is  found  that  ([2],  Part  2,  Sec.  VIII)  this  detection  probability  is  specifically 


PD=^1  +  ©| 


J*) 

^-0-' 

V2 


V 


l-aaf>7 


,  (Eq.  (69),  [2]  , 


with  the  associated  (conditional)  false  alarm  probability 


(5.7) 


GCp  =  ■ 


1-© 


of1  log  (K/ It) 
_2V2  V2o<*> 


U  =  p/q;  (Eq.  (71),  [2], 


(5.8) 


Here  is  the  detection  parameter  governing  performance,  given  by 
(*)  *  * 

C} '  =  varHo  gj  =  -2B) '  ,  cf.  (5.4a)  and  (5.4b)  in  the  coherent  cases) .  (5.9) 

The  quantities  K  and  (i  (=p/q)  are  respectively  the  threshold,  which  establishes  the  desired  false 
alarm  probability  dp  ,  and  the  a  priori  data  probability  ratio  ji  [14]. 


A.  Structure  of  the  Detection  Parameter,  ajp2: 

It  is  the  detection  parameter  we  wish  to  maximize,  or  at  least  make  as  large  as  practically 
possible,  by  choosing  suitable  input  signal  waveforms  so  that  the  combination  of  "output"  (i.e., 
received)  signal  structures  vis-k-vis  the  now  dominating  signal-dependent  noise  ('F(S),  (5.2)) 
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produced  by  reverberation,  cf.  Figure  5.1,  according  to  the  governing  geometry  of  propagation,  is 
made  as  large  as  possible.  For  the  coherent  regimes  considered  here,  (5.4b)  is  the  governing  form 
of  the  detection  parameter,  generally,  viz.: 


°i!=2L'2>J 


(5.10) 


where  s(nm)  is  the  real  normalized  [cf.  (5.4a)]  signal  received  at  the  m^  sensor,  at  time  tn,  while  a^ 
is  the  associated  (normalized)  amplitude,  cf.  (5.4f). 


Our  task  now  is  to  relate  the  real  received  signal  to  the  complex  injected  signal  j, 

in  order  to  establish  canonically  the  role  of  the  target  in  modifying  the  originally  transmitted  signal, 
cf.  Figure  5.1.  The  explicit  result  depends,  of  course,  on  the  target  model  chosen  as  well  as  on  the 
propagation  conditions.  From  Section  4,  we  have  the  desired  target  structure  factor  (TSF),  (4.7), 
(5.4f),  etc.  as  the  ratio  of  two  square-of-the-mean  signal  envelopes  in  these  coherent  detection 
cases: 


TSF :  S  |s<”’  (u)tf/p(iX 


tm,n  =t„-T0-e  +  Atm 


=  ;  l  =  tn  -T0  - e  ;  T0  =  |RoT|/c0  +|RoR|/c0  .  J 


(5.11) 


Here  T„  is  the  path  delay  between  source  and  receiving  array  (Ox  — » 0R) ,  and  the  delays  Axm  are 
specifically  given  by 

^.=(i.-U)  r./c„  ,  (5.11a) 

A  A 

where  i0,  ioRare  respectively  the  wavefront  normals  of  the  received  signal  and  beam  direction  (of 

the  maximum  lobe  of  the  adaptively  formed  beam  here),  [13].  The  result  for  s(nm\  the  (real) 
normalized  signal  at  the  m^1  sensor  at  time  tn,  is  found  to  be  (Sec.  9.1,  [1]) 


(5.12) 


(5.12a) 

In  this  way  we  have  identified  the  various  components  of  the  normalized  received  signal  with  the 
original  transmitted  signal,  so  that,  with  'F  =  *FA+G  +  *F(S)  of  (5.2)  and  (5.2a)  above,  we  see  that 
the  detection  parameter  (5.10)  now  becomes  explicitly  (albeit  still  in  canonical  form) 


4”’  =  Re{V2s„(tn  -  T„  exp{io>0(t0  -  T.  -  e„)}}  =  s„  ; 


A(m)2  A  2F(m) 
_aa _ —  rq-n 

2T  2¥ 


••a^=a0V«. 
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respective  structure  factors  of  the  various  reverberation  components  and  of  the  target  signal  (cf.  4.7) 
(all  of  which  depend  on  the  input  signal!),  and  the  quantity  Axm  appearing  in  tmn  ,  cf.  (5.11)  et 
seq.,  and  (5.13).  We  shall  refer  to  the  explicit  forms  of  these  structure  factors  in  B  ff. 

For  preformed  beams  we  formally  drop  the  dependence  on  (m)  in  the  above  (5.10)-(5.13), 
so  that  (5.13)  reduces  to  the  simpler  result 


CT(‘)  I  —  2li^N> 

a-TU  f 

To  Eol 

2F,(tn-T.-e.) 

°N-coh|  preformed  “  ^ 

{fs+fv+fb+>pA40/a?}  . 

2N 

(5.14) 


The  optimum  threshold  detection  algorithm  (5.3)  now  becomes 


with 


II 

gN(X)L  =  ^N-coh  +  l°g  H  *  X/(Xn)(0n>  »  X  =  {*„  =  X(0}  » 

n»l 

<$L  =  -2bL'L  =  Eq.  (5.14),  etc. 


(5.15) 

(5.15a) 


B.  Components  of  the  Channel  Model 

We  consider,  again  in  summary  form,  the  various  channel  components  that  we  need  for 
evaluating  the  detection  parameters  (5.13),  (5.14).  These  are  of  two  classes:  the  target  structure 
factor  (TSF),  for  both  adaptive  and  preformed  beams,  and  the  different  reverberation  structure 
factors  (RSFs),  all  of  which  are  distinguished  here  by  being  input  signal-dependent. 

From  (4.1),  extended  to  include  the  original  array  elements  (m),  we  write  generally  in  these 
coherent  detection  cases  for  the  TSF: 
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I.  TSF-adaptive  beamforming: 

'11{cxP{'ik"(2a"  'v<)(‘. +  At~  "  e.  -  *)}) 

d 

•(qL(4  tn  -T0  +  Axm  -  x))s0(x-e0)h 

(5.16) 

with  e0  to  maximize  ;  tn  =  Eqs.  (6.10a,b)  [1],  for  the  appropriate  Doppler  scenario;  k0  =  coo/co  . 

Here  G^  =  G^(f0)  for  these  narrowband,  far-field  cases  is  given  by  Eq.  (4.1a),  and  qL  is  the  local 
target  kernel,  (4.3).  For  bistatic  and  monostatic  regimes,  2a0,  (4.2a),  is  given  specifically  by 
(6.14a, b),  [1];  see  Figures  4.1  and  42.  Here  e0  is  chosen  to  maximize  the  signal  for  these 

A 

coherent  cases,  which  usually  means  picking  a  maximum  of  the  "carrier"  component  of  Re  ; 

A 

So.*  ,  of  course,  is  the  (complex)  envelope  of  the  input  signal,  as  discussed  in  Secs.  4,  8  [1]. 

Finally,  #4.  T  is  the  (complex)  beam  pattern  of  the  transmitter. 

The  relation  ((5.16)  for  the  TSF  simplifies  for  the  case  of  preformed  beams.  Thus,  we  have 
Eq.  (6.17a),  [1], 

II.  TSF -preformed  beams: 

*  0  formally 

^RT 

The  details  of  (5.17)  are  developed  in  Sec.  6,[1],  with  various  specific,  elementary  target  models 
presented  briefly  in  Sec.  7,  [1], 

In  a  similar  fashion  we  may  describe  the  Reverberation  Structure  Factors  (RSFs),  cf.  Sec. 

10.2  of  [1].  We  have 

HI.  Reverberation  Structure  Factors:  Fs  v  B  : 

We  have  normalized  our  results  with  a  mean-square  noise  ,  (2.3),  which,  however,  is  not 
stationary,  as  in  the  usual  treatments  [1],  [5],  v.  Thus, 

*  =  Kx (t, ,  t. ) + *A+a  +  V(tl ) ;  ¥(S)  =  Re  Kx (t. ,  t, )  S  Kx (t„ t, )  (5.18) 


f<K}  ->Fo-n  =Fo(tn  -T„;  eo,  ...) ;  Axn 

«4.t  " 


-in  (4.1).  (5.17) 
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from  (8.4),  [I],  where  X  =  reverberation  process  and  Kx  is  its  (complex)  covariance;  see  Eq. 

(5.20).  Since  our  signals  [Sec.  8.2,  [1],  and  Figure  5.1 ]  are  pulses,  of  finite  duration  T<(<T0),  it 
is  appropriate  in  the  first-order  treatment  here  to  select  the  path  delay  t,  =  T0  +T/2  =  Tc',  namely, 
the  midpoint  of  the  pulse  about  the  target  location  (at  Ol>  cf.  Figure  6J[lj).  Accordingly,  from  the 
fact  that  %  =  ^A^(FS  +  Fv  +  Fb)  =  KX(T',T')  now,  we  may  use  (3. 1 1)  to  represent  the 
reverberation  structure  factors,  generally,  by 

FR.s.v.,(To)  =  Re{Kx(T^,ro)/(A5)}  =  KIt/(Aj)  ,  (5.19) 

where  for  these  n.b.  input  signals  Kr(Tq,Tq)  is  given  by  (3.1 1),  with  t(=  t2)  — >  T'  and  with 

Fs  +  Fv  +  Fb  =  XrKii(T0',T:)/(aJ)  (5.20) 

to  be  used  in  (5.13),  (5.14).  The  first  term  of  (3.11)  embodies  the  (usually)  primary,  first-order  or 

single-scatter  contributions,  while  the  second  term  contains  all  higher  scatter  interactions.  The 

* 

signal-dependent  character  of  KR  is  at  once  evident,  where  here  is  described  specifically  in 
Section  8.2,[1]  for  the  present  applications.  As  we  noted  in  Section  3.2,  it  is  in  principle  possible 
to  calculate  some  of  the  scattering  functions,  a,  and  the  cospreading  functions,  dip.  The  critical 
role  of  geometry  appears  in  the  structure  of  qs,£i^,  as  expected.  Equation  (3.13)  is  an  explicit 
example  of  <js(t,v). 

Various  approaches  giving  explicit  results  for  Fr  are  listed  below,  with  some  comments: 

A 

Reverberation  Models,  KR 

FOM  Theory  [8]  <-  WSSUS  FR=S,V,B-+ 

Eqs.  (9.18),  (9.18a) 

Eqs.  (9.19),  (9.21):  "short  signals" 

Volume  (V),  surface  (S),  and  bottom  (B): 

Volume:  gauss,  nongauss  process 

Surface:  ocean  surface  (gauss);  non  gauss  at  small  angles,  N.B. 

Bottom:  nongauss;  ice  bottoms  also. 

These  scatter  channels  are  not  entirely  WSSUS  (wide  sense  stationary  uncorrelated  scatter)  channels. 


s.d.  J 

Remarks : 
(V) 
(S) 
(B) 


"  Classical "  [10] 
Eq.  (51):  Mx.<x>(t) 

Eq.  (50):  Mx(x) 
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C.  Optimization  by  Choice  of  Signal  Waveforms:  Remarks 

Finally,  an  important  question  here  is  the  choice  of  input  signal  waveforms,  to  optimize 
further,  or  at  least  increase,  the  detection  parameter  a) ' .  The  reason  for  this  is  that  increasing 
signal  level  alone  buys  little  improvement  in  performance  because  of  the  corresponding  increase  in 
the  reverberation,  which  is  signal-dependent.  Therefore,  one  must  seek  signal  waveforms  for  a 
given  signal  energy,  which  because  of  their  coherent  structure  vis-a-vis  the  reverberation,  can 

discriminate  against  the  latter,  in  effect  increasing  a} In  these  instances  the  customary  energy 
constraint  (used  against  an  additional  signal-independent  background  noise)  must  be  replaced  by 
one  that  preserves  the  area  under  the  ambiguity  function  of  the  input  signal,  or  requires  it  to  be  some 
specified  function  of  unity,  subject  to  bounds  on  v  and  x,  cf.  (3.9),  (3.10),  so  that  the  reverberation 

intensity  Kx(ti,tj)  =  'Fs+V+B  =  'P(S)  is  reduced  vis-^-vis  2^Json|  F0  in  (5.13),  (5.14),  for  example. 

Some  typical  signal  structures  for  this  purpose  are  described  in  Section  8.2  of  [1],  along  with  a 
sketch  of  a  formal  variational  procedure  for  maximizing  the  detection  parameter. 

6.  CONCLUDING  REMARKS  AND  NEXT  STEPS 

In  the  preceding  sections  we  have  provided  a  guide  to,  and  review  and  overview  of,  the 
analytic  methods  and  principal  results  of  the  author's  recent  approaches  [1]  to  the  treatment  of  the 
actual  threshold  detection  of  weak  targets  in  a  reverberation  dominated  ocean  environment  The 
aim  here  has  been  to  present  this  material  in  a  logical  sequence  and  in  sufficient  detail  so  that  the 
reader  can  proceed  from  this  work  to  implement  the  various  next  steps  (noted  in  part  below), 
which  are  needed  in  the  practical  application  of  the  analysis,  as  well  to  explore  and  use  the 
analytical  details  developed  in  the  original  study  [1].  In  this  sense  [1]  is  the  in-depth  "backup"  to 
this  paper. 

Here  we  have  summarized  our  results  for  the  following  elements,  all  of  which  play  then- 
significant  role  in  achieving  practical,  near-optimal  algorithms  and  performance  for  these  active 
detection  systems  in  typical  ocean  milieux  (cf.  Figure  5.1).  These  include 

(1)  Channel  models,  whose  component  processes  are  required  in  order  to  implement  the 

target  detection  algorithms  and  to  predict  performance  (Sections  2, 3); 

(2)  In  particular,  the  various  Structure  Factors  (cf.  Section  5.2)  needed  for 

A.  The  target  model(s)  -  which  here  are  spatially  linear  (cf.  Section  4); 

B .  Reverberation  -  General  results,  for  surface  (5),  volume  (V),  and  bottom  (B) 

(cf.  Section  3); 

C.  WSSUS  Channels  in  particular  (Section  3); 

D.  Multiple  scatter  contributions  (non-WSSUS  components),  Section  3; 
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The  analytic  models  for  these  are  specialized  to  account  for  the  principal  conditions  of  operation 
here,  namely: 

n.b.  =  "narrowband"  or  "short"  signals  ; 

N.B.  s  "  narrow -beams,"  usually  preformed ; 
f.f.  =  "far -field"  operation,  e.g.,  in  the  Fraunhofer  region ;  and 
s.d.  =  small  dopplers,  a  general  condition , 
with  independent  noise  sampling  and  zero-gradient  (Vc  =  0)  propagation  models. 

(3)  General  and  specific  (i.e.,  coherent)  threshold  detection  algorithms  (cf.  Section  5); 

(4)  Canonical  waveforms  and  noise  models,  both  ambient  and  reverberatory  (cf. 

Sections  2-4); 

(5)  Performance  measures,  both  canonical  and  specific  to  these  signal-dependent 
situations  (cf.  Section  5);  and 

(6)  Results  for  both  adaptive  and  preformed  beams. 

Some  Next  Steps: 

A  variety  of  next  steps  is  needed  (also  noted  in  [1].)  We  list  (not  necessarily  in  order  of 
importance)  the  following: 

(1)  Treatment  of  dependent  noise  samples,  primarily  spatial  [16]; 

(2)  (Vc  *  0) :  this  is  the  usual  situation  ([8],  parts  HI,  IV); 

(3)  Incoherent  detection  modes;  mixed  coherent  and  incoherent;  with  generalizations  of  the 
usual  ambiguity  function; 

(4)  Detailed  evaluation  of  Fs,  FR,S  V  B,  for  "short  signals"; 

(5)  Experimental  design  -  role  of  multiple  scattering  and  measurement  of  departure  from 
the  WSSUS  state  (cf.  Section  3); 

(6)  Determination  of  the  surface  structure  functions,  Fs,  in  detail; 

(7)  Adaptive  arrays,  "matched  field"  processing; 

(8)  Extension  of  the  general  n.b.  treatment  to  "broadband"  (f0/AF  <  2  or  3)  cases; 

(9)  Multiple  "ping"  integration; 

(10)  Comparisons  of  performance  for  different  signal  classes; 

(11)  Analytical  derivation  of  the  optimal  waveform  class  in  these  signal-dependent 
backgrounds;  and 

(12)  Estimation  ("classification")  problems  in  these  detection  environments  [16]. 

Some  additional  results,  and  next  steps,  are  also  cited  in  [1],  Sections  2,  2.2;  Section  3.G; 
Section  5;  and  Section  7.1. 
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